HOLOMORPHIC EXTENSION OF SOLUTIONS OF SEMILINEAR 

ELLIPTIC EQUATIONS 



MARCO CAPPIELLO AND FABIO NICOLA 

Abstract. We prove, for a wide class of semilinear elliptic differential and pseu- 
dodifferential equations in M. d , that the solutions which are sufficiently regular 
and have a certain decay at infinity extend to holomorphic functions in sectors 
of C d , improving earlier results where the extension was shown for a strip. More- 
over, exponential decay for such extended solutions is also proved. The results 
apply, in particular, to solitary wave solutions of many classical nonlinear evolu- 
tion equations. 



I. Introduction 

The main concern in this paper is the study of holomorphic extensions of the 
solutions of semilinear elliptic equations in Broadly speaking, we deal with 
equations of the form 

(1.1) Pu = F[u], 

where P is a linear elliptic differential, or even pseudodifferential, operator in M. d 
and F[u] is a nonlinearity, possibly involving the derivatives of u. For a wide class 
of equations of this type it is known that every solution u sufficiently regular and 
with a certain decay at infinity, actually is analytic on M. d and it extends to a 
holomorphic function in a strip of C d of the form 

(1.2) {z = x + iy eC d : \y\ < e}, 
for some e > 0, satisfying there the estimates 

(1.3) \u(x + iy)\ < Ce" cN 

for some C > 0, c > 0. A pioneering work on this subject was the paper by Kato and 
Masuda [23] . Later the problem of the holomorphic extension in a strip has been 
intensively studied in connection with the applications to solitary wave equations. 
In particular it was noticed in dimension 1 that several model equations like the 
Korteweg-de Vries equations and its generalizations, Schrodinger-type equations, 
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long wave-type equations admit solitary wave solutions which extend to meromor- 
phic functions with poles out of a strip of the form (jl.2p and having a decay of type 
(11. 3p . Among the main contributions in this sense we recall the papers by Bona 
and Li [8], [9], Grujic' and Kalisch [21], Bona, Grujic' and Kalisch [6], Bona and 
Weissler [10J. We also recall the paper [22] by Hayashi on well-posedness of the 
generalized KdV equation in Bergman-Szego spaces of holomorphic functions in a 
sector. Apart from its interest "per se" in the analysis of regularity properties of 
the solutions of differential equations, the study of complex singularities of solitary 
waves could give information also on the onset of real blowup; we refer to Bona 
and Weissler [10J for a fascinating discussion in this connection. 

Recently, the properties described above have been proved for some general 
classes of semilinear elliptic equations in any dimension, even with variable coeffi- 
cients; see for example Biagioni and Gramchev [5], Gramchev [20] and Gramchev, 
Cappiello and Rodino [HI [121 EI21 US]- The results in these papers have been 
stated and proved in terms of estimates in the Gelfand-Shilov spaces of type S, cf. 
[19]. We refer to Nicola and Rodino [271 Chapter 6] for a self-contained account of 
these results. 

Nevertheless in the above mentioned papers some relevant issues remained un- 
explored. The first one is the identification of a maximal holomorphic extension. 
In other words, the problem is to understand what is the biggest complex domain 
on which a holomorphic extension is possible. The second one is related to a dual 
aspect and concerns the identification of the maximal domain on which the decay 
properties on M. d of solutions remain valid. 

To be precise, let us introduce a class of operators which the above results and 
those in the present paper apply to. First, one can consider differential operators 
with polynomial coefficients 

(1.4) P= Yl c ^D a , 

\a\<m, \p\<n 

m > 1, n > 0, c a/ 3 G C, with symbol 

(1.5) p(x,s)= Yl c ^ x ^ a 

\a\<m, \P\<n 

of G-elliptic type, namely satisfying the following global version of ellipticity: 

(1.6) Yl c ^ x ^ a ^ °> ^ R '.^°. 

\a\=m, |/3|<n 

(1.7) Yl c ^ x ^ a + °> ^°.^ R - 

\a\<m, \P\=n 
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(1.8) Yl c ^ x ^ a ^ °> x ^ °> ^ °- 

|«|=m, |/3|=n 

If P has constant coefficients (that is n = 0) then (11.61) . (11.71) . (ll.Sp are satisfied if 
and only if P is elliptic and its symbol p(£) satisfies p(£) ^ in M. d . As a relevant 
model one can consider the operator P = —A + A, A 6 C, A ^ R_ U {0}. 

More generally we deal with pseudodifferential operators. Namely, given real 
numbers m > 0, n > we consider symbols p{x,^) satisfying the following estimates 

(1.9) \dgd%p(x,Z)\ < C H+m+1 a\(3\{0 m ' H (x) n - m 

for every (x, £) G M. 2d ,a,(3 G N d and for some positive constant C independent 
of a, (3 (we denote as usual (x) = (1 + Is) 2 ) 1 / 2 ). Notice that every polynomial 
p(x,£) as in (II. 5p satisfies such estimates. For this general class the condition of 
G-ellipticity can be stated by requiring that 

(1.10) \p(x,0\>c(O m (x) n , \x\ + \£\>C, 

for some positive constants c, C. It is known that a symbol of the form (11.51) satisfies 
(II.IOP if and only if it satisfies the three conditions (11.61) . (II. 7p . (11.81) simultaneously. 

Finally, concerning the nonlinear term F[u] in (11.11) . we assume here that it is of 
the form 

i 

(i-ii) **[«] = £ E F h,i >P u..., P ^ h Il dPk ^ 

hi pi,...,pi k=l 

where h G N d , with < \h\ < max{n-l, 0}, p k G N d with < \p k \ < max{m-l, 0}, 
/ G N, / > 2 and Fh,i, Pl ,..., Pl G C (the above sum being finite). Moreover, we allow 
some of the factors in (11.111) to be replaced by their complex conjugates. 
As a simple example, consider the equation 

(1.12) - Au + u = |u|^V 

with / G N, / > 2 odd, which arises e.g. when looking for standing wave solutions 
for the Klein-Gordon or Schrodinger equation, as well as travelling wave solutions 
for the Klein-Gordon equation (cf. Berestycki and Lions \A\). The existence of 
solutions in /7 1 (M a! ) and their exponential decay was studied in [1], whereas the 
possibility of extending them holomorphically on a strip has been recently shown 
in [13] (incidentally, the exponential decay generally drops for elliptic equations 
which are not globally elliptic, such as —Am = |w| z-1 u). 

The above more general class of operators includes, in dimension d = 1, the 
solitary wave counterpart of several evolution equations of Korteweg-de Vries type, 
as well as of long-wave type; see [9] and below. This also motivated the statements 
of the results for pseudodifferential operators, or at least for Fourier multipliers 
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(e.g. p(£) = £Coth£ + A, A > —1, for the intermediate-long-wave equation; see 

Now, it is known from [9] and [T3J Theorem 7.3] that (for classes of nonlinearities 
that intersect the above one), if P is G-elliptic, then all the solutions u of fll.l[) with 
u E H s (R d ), s>| + m- lforn>0 (respectively (x) £o u E H s (R d ), s > | + m — 1 
for n = 0) actually decay at infinity like e _c ' x ',c > 0, and extend holomorphically 
on a strip of the form ( 11.2H . 

Here we shall improve this result by showing that the holomorphic extension and 
the exponential decay actually hold in a sector of the complex domain. Namely, we 
have the following result, which seems new even for the simplest equation (I1.12p . 

Theorem 1.1. Let P be a pseudodifferential operator with symbol p satisfying 
(11. 9p . m > 0, n > 0, and assume that p is G- elliptic, that is (ll.lOp is satisfied. 
Let F[u] be of the form (11.111) (possibly with some factors in the product replaced 
by their conjugates). Assume moreover that u E iP(IR d ) ; s > | + max{|p&|} ; is 
a solution of (II. ip . In the case n = assume further (x) £0 u E L 2 (JBL d ), for some 
eo > 0. Then u extends to a holomorphic function in the sector of C d 

(1.13) C E = {z = x + iy E C d : \y\ < e(l + \x\)}, 

for some e > 0, satisfying there the estimates (II. 3p for some constants C > 0, c > 0. 

The width of the sector depends in general on the solution u considered; in fact, 
for autonomous equations the width must depend on u, because one can exploit the 
invariance with respect to translations to move the complex singularities parallel 
to R d . 

The shape of the domain of holomorphic extension as a sector is sharp, in the 
sense that, even in dimension d — 1, for any angle 6 ^ 0, tc we can find G-elliptic 
equations with constant coefficients admitting Schwartz solutions whose meromor- 
phic extensions have a sequence of poles along the ray arg z = 9. We refer to 
Section below for details on this point and also for remarks on the a priori regu- 
larity assumptions in Theorem 11.11 

The linear case (F[u] = 0) deserves a separate discussion. Indeed, the analysis 
of the linear equation Pu = is important for the holomorphic extension and the 
decay of eigenfunctions of G-elliptic operators and their powers; see Maniccia and 
Panarese [25] and Schrohe [31]. Moreover in the linear case it is possible to relax 
the assumptions on the regularity and decay of u, admitting solutions with a priori 
algebraic growth and in that case the width of the sector is independent on the 
particular solution considered; see Theorem 15.11 below. 

Finally we present an application of the above result to solitary waves. Following 
P], we consider the following class of Korteweg-de Vries-type equation 



(1.14) 



vt + v x + F[v] x - (Mv) x = 0, (t, x)ERxR, 
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and long-wave-type equations 

(1.15) vt + v x + F[v] x + (Mv) t = 0, (t,x)eRxR, 

where M = p(D) is a Fourier multiplier, F[v] is a polynomial with real coefficients, 
F(0) = F'(0) = 0, and subscripts denote derivatives. We look for solutions v(t,x) 
of solitary wave type, i.e. v(t, x) = u(x — Vt), for some function u of one variable 
and some constant velocity V. We have the following result. 

Theorem 1.2. Let satisfy the analytic symbol estimates of order m e 1, 
namely 

(1.16) \d a p{0\ < A a+1 a\(O m ~ a , (61, «GN 
for some constant A > 0, as well as the lower bounds 

(1.17) P (0 >0,(el, p(0 > C-^r, |£| > C, 

/or some constant C > 0. Suppose moreover m > 1. 

Lei f(t,x) = m(x — Vt) fre a u>eaA; solution of (11. 14)) or (11.151) . V > 1, 
w G L°°(]R), lim^-s.oo m(x) = 0. T/ien w extends to a holomorphic function u(x + iy) 
in the sector 

(1.18) {z = i + i|/6C: 1 2/ 1 < e(l + |x|)} 

/or some e > 0, satisfying there the estimates (jl.3p /or some constants C > 0, c > 0. 

Notice that the estimates (I1.16P are satisfied by any polynomial p(£) of degree 
m. More generally, the condition (11.161) is equivalent to saying that p(£) extends to 
a holomorphic function + iff) in a sector of the type (11.181) . and satisfies there 
the bounds |p(£ + ^)| < C'(^) m (see Proposition 15.21 below). This remark makes 
(11. 1 6 p very easy to check in concrete situations, where typically p(£) is expressed 
in terms of elementary functions. 

We also observe that, under the hypotheses of Theorem 11.21 we already know 
from [5] that u extends to a holomorphic function on a strip and displays there 
an exponential decay of type (11.31) . hence Theorem 11.21 can be regarded as an 
improvement of this result. For the existence of solitary waves for these equations 
we refer to the detailed analysis in Albert, Bona and Saut [lj, Amick and Toland 
[2], Benjamin, Bona and Bose [3], Weinstein [33] . 

Finally we mention that similar extensions results should hopefully be valid for 
other classes of non-linear elliptic equations, e.g. with linear part P = —A + 
| x | 2 (cf. [llj). Similarly, even non-elliptic hypoelliptic operators and dispersion 
managed solitons could share similar properties, possibly with the sector replaced 
by a smaller domain (but larger than a strip). We plan to investigate these issues 
in a subsequent paper. 

The paper is organized as follows. Section [2] collects notation and some prelim- 
inary results about G-pseudodifferential operators (composition, boundedness on 
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Sobolev spaces, parametrices, etc.). In Section [3] we introduce a suitable space of 
analytic functions on ~R d , which admit a holomorphic extension to sectors in C d , 
and we prove some relevant properties used in the sequel. In Section @] we prove 
Theorem 11.11 The proof is based on the application of an iterative Picard scheme in 
the space of analytic functions defined in Section [31 In Section [5] we prove Theorem 
11.21 and treat in detail the linear case F[u] = 0, see Theorem 15.11 We moreover 
show some other examples and counterexamples and test on them the sharpness of 
our results. 



2. Notation and preliminary results 

2.1. Factorial and binomial coefficients. We use the usual multi- index nota- 
tion for factorial and binomial coefficients. Hence, for a = (cti, . . . ,ctd) £ N d we 
set a! = a%\ . . . and for /3, a £ N d , f3 < a, we set Qj) = p^-py ■ 
The following inequality is standard and used often in the sequel: 

(2.1) ( a ] < 2H 







Also, we recall the identity 



\a'\=3 

ClJ <-OL 

rd 



which follows from ULii 1 + = i 1 + *) H > and 

gives in particular 



(2.2) Q*Q- 

The last estimate implies in turn, by induction, 
(2-3) TT""~~T7 — TTTi^TFTT ' « = *i + ■ ■ ■ + * 



as well as 

{2A) t^W^R^' 

Finally we recall the so-called inverse Leibniz' formula: 

(2.5) x p d a u{x)= 7WZ-§ ( a )d a ^(x^u(x)). 
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2.2. G-Pseudo-differential operators. Pseudo-differential operators are formally 
represented as integral operators of the type 

(2.6) p(x, D)u{x) = (2vr)- d / e^p(x, £)£(£) 
where 

u(Q = / e~ lxi u(x) dx 

denotes the Fourier transform of u and p(x, £) is the so-called symbol of p(x,D). 
According to the symbol spaces which p belongs to, one can consider u in several 
classes of functions or distributions and symbolic calculi and boundedness results 
on Sobolev spaces are available. We briefly recall this for the class of the so-called 
G-pseudodifferential operators (also named SG or scattering pseudodifferential op- 
erators in the literature). They are defined by the formula (12.61) . where p(x, £) 
satisfies, for some m,n G M, the following estimates: for every a, (3 G N d there 
exists a constant C a ^ > such that 

(2.7) \d%dPp(x,Z)\ < C a ^(x) n -W(O m - H 

for every x, £ G M d . The space of functions satisfying these estimates is denoted by 
G m,n (IR d ), whereas we set OPG m,n (IR d ) for the corresponding operators. We endow 
G rn ' n (M d ) with the topology defined by the seminorms 

||p||i?= sup sup {\d^p(x,0\(x)~ n+m (0~ m+H }, NeN. 

|et|+|/8|<JV (x,C)eM 2d 

As a prototype one can take P = —A + A, A G C, of order m = 2, n = 0. 
More generally, the case of Fourier multipliers, where the symbol p(£) depends 
only on £ (hence n = 0) is of great interest, mostly for applications to solitary 
waves and ground state equations, see [13J. As another example, we have x^d 01 G 
OPGl a l^l(R d ). 

The classes OPG m,n (M. d ) were introduced in [28] and studied in detail in |16j . 
[T7] . [26], [30]. They are in fact a particular case of the general Hormander's classes, 
see [231 Chapter XVIII], and turn out to be very convenient for a series of problems 
involving global aspects of partial differential equations in M. d . 

We now summarize some properties which will be useful for us later on; beside 
the above mentioned papers, we refer to [23, Chapter 3] for a detailed and self- 
contained presentation. 

First, if p G G m ' n (M. d ) then p(x,D) defines a continuous map 5(IR a! ) — > S(M d ) 
which extends to a continuous map <S'(M ) — > S'(WL d ). The composition of two 
such operators is therefore well defined in 5(IR d ) and in 5'(M d ); more precisely, 
if pi G G mi ' ni (R d ) and p 2 G G m2 > n2 (R d ), then p t (x, D)p 2 (x, D) = p 3 (x,D) with 
p 3 G G mi+m2 ' ni+n2 (R d ) and the map (pi,p 2 ) Ps is continuous G mi ' ni (R d ) x 
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A symbol p G G m,n (R d ) (and the corresponding operator) is called G- elliptic if 
for some constants C, c > 0, it satisfies 

(2.8) \p(x,0\>c(x) n (O m , M + KI>c. 

For example, P = —A + A is G-elliptic if and only if A ^ R_ U {0}. More generally, 
for an operator with polynomial coefficients as in (jl.4p . G-ellipticity is equivalent 
to (Ojl . ffLTj) . (Oj) to hold simultaneously 

G-ellipticity guaranties the existence of a parametrix E G OPG~ m, ~ n (R d ) of 
P = p(x, D), namely EP = I + R and PE — I + R', where R, R' are (globally) 
regularizing pseudodifferential operators, i.e. with Schwartz symbols. Hence R and 
R' are continuous maps S'(IR d ) — >■ S(IR d ). 

It is useful to consider the action of such operators on the standard Sobolev 
spaces 

H s (R d ) = {ue S>(R d ) : ||«||. := Q |^(0| 2 (1 + \^) s < oo}, 

and on the weighted versions 

H Sl ' S2 (R d ) = {ue S'(R d ) : ||m|| Si , S2 := ||(x) S2 w|| Sl < oo}. 
Indeed, if p G G m ' n (R d ) then 

(2.9) p(x, D) : H Sl ' S2 {R d ) -»■ H Sl -' m ' S2 - n {R d ) 
continuously, and 

\\P(X, D)\\ B ( H s 1 ,s 2 ( R d^ H s 1 - mi ,s 2 -m 2 ( R dy ) < C\\p\\pf 

for suitable C > 0, iV G N depending only on s\, S2, m 1 ,m2 and on the dimension d. 
In particular, for s 2 = we see that, if n < then p(x, D) : H s (R d ) — > H s ~ m (R d ) 
continuously for every s6M. 

Remark 2.1. The complex interpolation for the spaces H Sl,S2 (R d ) works as one 
expects, i.e. for s±, S2, t±, t% G R, < 8 < 1, 

(2.10) [H Sl,tl (R d ) , H S2,t2 (R d )]e = H s ^(R d ), s = (l-6)si+6s 2 , t = {l-9)t x +9t 2 , 

see for example [Hl[32]. The property (I2.10p will be useful in the sequel in view of 
the following consequence: suppose that u G H s (R d ) and (x) £0 u G L 2 (R d ) for some 
So > 0. Then for every s' < s there exists e > such that (x) £ u G H s (R d ). 

Finally we point out for further reference the following formulas, which can be 
verified by a direct computation: for a, (3 G N d , 

(2.11) x p Pu = V(-1) M r J (Djp)(x, D)(x^u), 



7</3 V U 



HOLOMORPHIC EXTENSION FOR ELLIPTIC EQUATIONS 



(2.12) d a Pu = (fj (d S x p)(x, D)d a ~ s u. 

3. A SPACE OF ANALYTIC FUNCTIONS 

We introduce here a space of analytic functions in M d , already considered in [12] 
(and denoted there by 5*^*(lR d )), which is tailored to the problem of the holomorphic 
extension to subsets of C d of the form (11.131) . as shown by the subsequent Theorem 

Definition 3.1. We denote by *4 sec t(R d ) the space of all functions f E C°°(lR d ) 
satisfying the following condition: there exists a constant C > such that 

(3.1) \x^d a f(x)\ < C H+l/3|+1 max{H, \/3\}\, for all a, (3 E N d . 

Theorem 3.2. Let f E *4. serf (IR d ). Then f extends to a holomorphic function 
f(x + iy) in the sector of C d 

(3.2) C £ = {z = x + iy EC d : \y\ < e(l + |x|)} 
for some e > 0, satisfying there the estimates 

(3.3) \f(x + iy)\<Ce-M, 
for some constants C > 0, c > 0. 

Proof. First we show the estimates 

(3.4) \x p d a f(x)\ < C |a|+1 |a|!e- c|x| , for \(5\ < \a\. 

Indeed, since \x\ n < & n ^]| 7 | =n |^ 7 | for a constant k > depending only on the 
dimension d, by (13. 1 p we have (assuming C > 1 in (13. ip ) 

e c ^\x^d a f{x)\ = f2^Y^\^d a f(x)\ 

n=0 n ' 

oo 1 



n=0 \l\=n 



7 



oo 



(l«l 


+ l7l)! 




ItI 


i 



n=0 \~t\=n 

Since the number of multi- indices 7 satisfying I7I = n does not exceed 2 d+n_1 , an 
application of (I2.ip gives 03. 4p for a new constant C, if c is small enough. 
Now, 03. 4p and the estimate |a|! < <i' a 'a! give 

(3.5) \d a f{x)\ < C |Q|+1 a!(x)- |a| e- c|a;| , 



10 



MARCO CAPPIELLO AND FABIO NICOLA 



for a new constant C > 0. This shows that the power series 

d a f(x) 



(3.6) £ 



z - x 



a l 



for any fixed x G W 1 converges in a polydisc in C d defined by \zk — Xk\ < 
1 < k < d. The union of such polydiscs, when x varies in M d , cover a subset 
C £ C C d of the type (13.21) . for some e > 0. Since on the intersection (when not 
empty) of two such polydiscs these extensions agree {R d C C d is totally real), f\x) 
extends to a holomorphic function on C £ . For z G C £ , using the representation (13.61) 
as a power series with x = Rez and (13.51) . we also get the desired estimate (13.31) 
for a new constant C > 0. □ 

In the sequel we will use the following characterization of the space ^4 sec j(M^) in 
terms of i? s -based norms. 
Set, for / G S'(R d ), 



Proposition 3.3. Let f G A sec t(R d )- Then for every s > there exists e > such 
thatS s J[f}<oo. _ ' 

In the opposite direction, if for some s > there exists e > such that S^[f] < 
oo then f G A sec t( 



Proof. Assume / G A sec t0& )• It suffices to argue when s = k is integer. Then 

H^/iu < c £ ||^(x^./-)iu 2 . 

| 7 |<fc 

Now, if M G N satisfies M > d/iwe have 

(3.8) \\d^{x^d a f)\\ L 2 <C n \\(\ + \x\ 2 ) M &{xP& x f)\\ Lao . 

On the other hand we have, by Leibniz' formula and (13.11) . 

\\(i + \x\ 2 ) M d^(x^d a f)\\ Lao < (^jg^^ + l^^^^fh- 

< C 7 |/3| |7| C |Q|+I/31 max{2M + |a| + |7|}!. 

Since max{2M + \/3\, \a\ + | 7 |} < max{|/3|, |a|} + 2M + | 7 | and \/3\^ < C\ P \ by 
f l2~Tj) we get 

||x^/|| fe <Ci a|+l/3|+1 max{|«|,|/3|}! 
for some constant > 0. Hence S'^ff/] < oo if £ < . 
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In the opposite direction, we may take s = 0; hence assume S^[f] < oo for some 
e > 0. Then (13.11) holds with the L°° norm replaced by the L 2 norm. If M is an 
integer, M > d/2, we have 

\\x p d a f\\ L oo <C W{x (3 d a f)\\ L 2. 
h\<M 

Hence an application of Leibniz' formula and the same arguments as above show 
that / G A sect (R d ). U 



4. Proof of the main result (Theorem 11.11 ) 

In this section we prove Theorem 11.11 In fact we shall state and prove this result 
for the more general non-homogeneous equation 

(4.1) Pu = f + F[u], 

where P and F[u] satisfy the assumptions of Theorem 11.11 and / is a function in the 
space A se ct(R d ) defined in Section 3. Moreover we shall restate our result in terms 
of estimates in A sec t(M. d ). Namely, in view of Theorem 13. 2\ it will be sufficient to 
prove the following theorem. 

Theorem 4.1. Let P be a pseudodifferential operator with symbol p satisfying 
(II. 9j) . m > 0,?2 > and assume that p is G -elliptic, that is fll.lOf) is satisfied. Let 
F[u] be of the form (11.1 ip (possibly with some factors in the product replaced by 
their conjugates) and f G A sec t(R d ). Assume moreover that u G H s (W i ), s > | + 
max{|pfc|}, is a solution of (14. ip . In the case n = assume further (x) £ °u G L 2 (R d ) 
for some Eq > 0. Then u G A sec t(R d ). 

In fact we always assume that F[u] has the form in ( 11.111) . and we leave to the 
reader the easy changes when some factors of the product in (II. lip are replaced by 
their conjugates. 

We start by showing that, under the assumptions of Theorem 14. 1| u is in fact a 
Schwartz function. 

Lemma 4.2. Let P, m, n, F[u], f be as in Theorem \4-l\ Let u be a solution 
of (14. ip satisfying (x) £ °u G H s (R. d ) for some s > ~ + max{|p fc |}, £o > 0. Then, 
when n > 0, we have (x) £o+r72 (D) ai u G H s (M> d ) for every G\ < min{m, 1} ; and 
o~2 < min{n, 1}. If n = and we assume in addition Eq > 0, then we have 
(x) £o+a2 (D} ai u G H s (R d ) for every <j\ < min{m, 1} and cr 2 < min{e , 1}- 

Proof. We first consider the case n > 0. Let E G OPG~ m, ~ n (M> d ) be a parametrix 

for P; hence R:=EP -I G OPC^ 1 '- 1 ^). We have from (Ol) 

(4.2) 

(x) £0+a2 (D) ai u = (x) £0+a2 {D) ai Ef - (x) £0+a2 (D) ai Ru + (x) £0+a ' 2 (D) ai EF[u]. 
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Since a x < 1 and a 2 < 1, the operator (x) a2 (D) ai R G OPG°'°(IR d ) is bounded 
on H s ' £ °(R d ); cf. ( 12. 9p . Taking also into account the assumptions on / it follows 
therefore that the if s -norm of the first two terms in the right-hand side of ( 14. 2 j) is 
finite. Concerning the last term, observe that, by the assumptions on 01,02 and 
h, the operator (x) a2 (D) ai o E o x h G OPG- m+<Jl '- n+CT2+|/l| (M d ) belongs in fact to 
QP£-ma X {m-i,o},0(]g>d^ ^ g a conse quence, since M := max{|p fc |} < max{m-l,0}, 
it is bounded H s ~ M ' £0 (M. d ) — > H s ' £0 (R d ). Hence by Schauder's estimates (recall 
that s > d/2 + M) we have 

1 

\\(x) £0+ ^(D)^EF[u}\\ s <C s J2 E ll^} £o n^ w il^ 

Z pi,...,p ; fc=l 

M II s-M 

/ pi,...,pi k=l 
I 

^"E E nii^^ii^^^'Eii^ii^ 00 - 

i pi,...,pi k=l I 

We treat now the case n = 0, hence /i = in the nonlinearity ( II. lip . We consider 
again (@~2]). For the terms (x) £o+a2 (D) ai E f and (x) £o+(72 (D) ai Ru we argue as 
before. For the nonlinear term observe that, since M := max{|/?fc|} < maxjm — 
1,0}, we have that (D) ai E G OPG- m ^ m ~ 1 '°^°(R d ) is bounded H s - M ' l£0 (R d ) ->• 
H s ' l£ °(R d ), for every I; see (I2.9p . Hence for a 2 < e < e (l — 1) we get 

1 

\\(xr o+ ^(D)^EF[u]\\ s < C S ]T W(x) lE0 (DrEl[d^u\\ s 

l pi,.--,pi k=l 
I 

U II s-M 

I pi,...,p; fc=l 
/ 

^ ^e e n ii(^) £o ^^h s -m 

l pi,...,pi k=l 

< c':'J2\\(xy°u\\ l s <oc, 

I 

where we applied again Schauder's estimate and, in the last inequality, ( 12. 9 p to d Pk . 

□ 

Let us observe that, when n > 0, an iterated application of Lemma 14.21 shows 
that, under the assumptions of Theorem 14.11 (x) T2 (D) Tl u G H s (M. d ) for every T\ > 
0, r 2 > 0, that is u G 5(lR d ). The same is true when n = 0, because the assumptions 
u G H s (R d ), s > d/2 + max{|p fc |}, and (x) £ °u G L 2 (R d ), e > 0, imply that for 
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new values of s and Eo as above we have {x) £0 u £ H s (M. d ) (see Remark 12. ip . and 
Lemma 14.21 still allows us to upgrade regularity and decay. 
In particular, the sum S^fw] is finite for every iV £ N. 

In order to prove Theorem 14.11 it suffices to verify that S^fw] < oo, in view 
of Proposition 13. 31 This will be achieved by an iteration argument involving the 
partial sum of the series in (13.71) . that is 

(43) sm ~-^^^ m 

We shall treat separately the cases m > 1 and < m < 1, since the study of the 
nonlinearity requires different arguments. 

4.1. Proof of Theorem I4.lt the case m > 1. We need several estimates to 
which we address now. 

Proposition 4.3. Let R £ OP(j _1 ' _1 (R d ). Then for every s £ R there exists a 
constant C s > such that, for every e < 1, 

Proof. We first estimate the terms with a = 0, hence ^ 0. Let j £ {1, . . . , d} 
such that (3j ^ 0. Since R o Xj £ OPG -1 '°(R d ) is bounded on H s (R d ) we haveS 

j^\\R(x^u)\\ s <C s e £ -^\\x^u\\ s . 

Similarly one argues if = 0, a ^ 0. If finally a ^ 0, /3 ^ 0, hence for some 
j, k £ {1, ... , d}, we have o^- 7^ 0, 7^ 0, we write 

and use the fact that Rdj o x k £ OPG°'°(R d ) is bounded on H s {R d ). We get 



c-lal+l^l r M+|/3|-2 



||i?(^9 Q it)|| s < CV 2 ri 1 i Q n, ll^~ efc ^~ ej 



max{|a|, |/3|}! max{|a|, |/3|} 



£ |a|+|/3|-2 



+ ° s£ max{|a|-l, |/3|-1}! 11 l|s ' 

(we understand that the second term in the right-hand side is omitted if f3j = 0). 
These estimates give at once the desired result if e < 1. □ 



We denote by ej the jth vector of the standard basis of K d . 
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Proposition 4.4. Let P = p(x,D) be a pseudodifferential with symbol p(x,£) 
satisfying the estimates (|r5|) , with m>0,n>0. Let E e OPG-" 1 - n (R d ). Then 
for every seR there exists a constant C s > such that, for every e small enough 
and u G 5(IR d ) ; we have 

_. F \«\+\P\ 

<«> £ max{ | aM/i | }! w f '-^Mi.<^y- 1 M. 

0<|a| + |/3|<Af 11 l ' l ^ IJ 

Proof. We have 

[p, = [P, a^ff* + [P, d a ] . 
Hence, using (l2~TTj) . f[27L2]) . we get 

(4.5) [P,:/<9> = ^ (-l) l7o|+1 ( ^ J(i5j p)(^^)(^ 70 ^) 



0^<Ka 



Given ft, 8, let 5 be a multi-index of maximal length among those satisfying \8\ < \8\ 
and 8 < ft (hence \8\ = \8\ unless (3 — 8 = 0). Writing x 13 = x 6 x^~ s in the last term 
of (14. 5p and using again (12. lip we get 
(4.6) 

[P,x p d a )u = J2 (-l) l7o|+1 ^ ~ J (^x^Dfdlp^x.D^-^d^u). 

S^QL iQ<ft — 8 

(5,7o)^(0,0) 

We now look at the operator x P- s -~yoga-5 _ Qj ven ^ a ^ ^ j Q k e a multi-index 
of maximal length among those satisfying |7 | < |7o| and 70 < a — 8. We write, by 
the inverse Leibniz formula (I2.5p . 

(4 7) x^~^~ 10 d a ~ 6 = x^~ s ~ 10 <9 a_<5 ~ 70 = <9 70 o x^ 6 ^ d a ~ s ~^ !o 

+ (~l) l7l| (/3 -8- To)! /7o\ g70 - 7l q x ^5-y -y ld a-S-^o 

»JzU^ W -~ s ~io -71)! W 



0^71</3-*-70 
71<7Q 



We now look at the operator s 70 71 <9 a 5 70 . We denote by 71 a multi- index of 
maximal length among those satisfying |7i| < |7i|, 71 < a — 8 — % . Again by the 
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inverse Leibniz formula we have 

(4 8) x P-&-~IQ-~1i Qa-fi-lO — x P-&-lti-li Qli QaS-lO-ll 

— Qli a; /3-<5-7o-7i ( g a -'5-70-7i 

\ r ( — l)l 7 — £ — 7q — 7i)! / 7i\ rp l _^ 2 ^ ^ ) g_ ( 5_ 70 _ 7l _ 72 g a _g„^ _^ 1 

„ <rf (/3-5-70-71-72)! V7 2 ; 

0^72</3- 5 — TO — n 111/ 
72<7l 

Continuing in this way and substituting all in (14 .6p we get 

h 

[P,xPd«]u = J2H E E • • E ^-70,7! 7i 

5< Q j=0 70<£~'5 0#7l</9-«-70 0^t j -<jS-J-7o-...— Yj_i 
(<5,70)^(0,0) 71<70 7j<7j_l 

where 7, is defined inductively as a multi-index of maximal length among those 
satisfying |tj| < |tj| and jj < a — 5 — % — ... — jj-i, 



lr , a\(J3-8)\ yr 

' a,/3A7 °' 7l '-' 7jl (a - 5)\5\ l0 KP - 5 - 70 - ... - 7,)! l\ V 7, 

(4.9) < H!|fi-£|! 2 l^+-^- 1 l, 

|a — <5|!5!7o!|/3 — 5 — 70 — ... — 

cf. (E3D and (ED), and 

(4.10) p QlM70 , 71 ,..., 7j (x,0 = ^(/v^)(,-.oC : " ~ • ' ' ■ J > 0, 

(if j = in (14. 9 p we mean that there are not the binomial factors, nor the power 
of 2). Observe that, since we have 7^ 7^ for every j > 1, this procedure in fact 
stops after a finite number of steps. 

Now we observe that, by (II. 9p . (12. ip . and Leibniz' formula, for every 6, a e N d 
we have 

(4.11) |^9> aiW , 70 , 7li ..., 7j (x,0l < ^"l + l 5 l +1 7o!5!(x)"- |ff| (O mH91 , 

for some constant C depending only on 9 and a. In fact \8\ < \8\, I70— 71+71 — . • • — 
lj + lj\ < \lo\ < |7o|j an d the powers of \8\ and |7o| which arise can be estimated 
by CItoI+H+i f or some c > 0. 

We now use these last bounds to estimate E op a ^{ i7l)i71] ... i7J (ir, D). To this end, 
observe that this operator belongs to OPG '°(lR £i ), and therefore its norm as a 
bounded operator on H s (M> d ) is estimated by a seminorm of its symbol in G 0,0 (R d ), 
depending only on s and d. Such a seminorm is in turn estimated by the product 



16 MARCO CAPPIELLO AND FABIO NICOLA 

of a seminorm of the symbol of E in G- m ^ n (R d ) and a seminorm of Pa,/3A7oi7i)--->7.j 
in G m ' n (M> d ), again depending only on s, d. Hence we see from (14.111) that 

(4.12) pop«^ATD,-n....,7 J (^.^)llB( J ff-(M-)) < ^l + l 5 l +1 7o !5!. 

Let now \f}\ > |a|. Then |5| = Using moreover the estimate 

\a\\\P-~5\\ _ \a\m-\W 
|/?|!|a-*|! |/3|!(|a|-|5|)! 

together with (14. 9[) and (I4.12p . we obtain 



< 1, 



(4.13) 



\C a , PAl0 ^... n] \\\Eop^^ 



u 



s 



< C s {C s e)\ & \ 



+|70+-+7il+]70+...+7j-l| 



|a|+| / 8|-|5|-|7o+...+7i|-|7o+-+7j-i| 



|/3- 8- 70 - ... - 7j|! 

X ||^-5-70-...-7^a-a-70-...-7 Jlt | 



Similarly, if \a\ > \(3\ we have 1%! = |7fc|, < k < j, and 

(414) |^-g|!|a- < y-7o---7 J -|! <1 

|a-5|!|/3-5- 7 o-...-7 J |! " ' 

(recall that if \5\ < \5\ then — 6 — 70 = . . . = jj — 7q = . . . = 7j = 0). 
By flUD, (gJ2D (|Q3jl . we get in this case 



(4.15) 

1/31 



e [«M0l 



a ! 



< (7 s ((7 s £)l 5 W7o+---+7jl+l7o+---+7j-il 



£ |a|+| / 9|-|5|-|7o+...+7 3 -|-|70+-+7j-l| 
|a - 5 -70 - . . . - 7j|! 
x 1 1 x p ~' s ~' 10 ~ ■ ■ ■ -~iW a - 5 -* 10 - ■ ■ ■ u I 



Since, if > \a\, we have 

max{|/3 - 6 - 70 - . . . - 7j|, |ar - 5 - 70 - . . . - 7,-|} = \/3 - 5 - j - . . . - 
whereas if |a| > \(3\ it turns out 

max{|/3 - 6 -70 - ■ ■ ■ -7j-|, \a- S - . . . = \a - S -70 - . . . - Til, 
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we deduce from ( 14. 13ft and ( 14. 15ft that, if e is small enough, 

^ £ \ a \ + \0\ 

M + |/3|<7V 11 h|M|J 

^ max{|«|,|/3|}!" ^ V 

\a\ + \P\<N-l 11 l,l ^ IJ J=0 <5 7i #0,..., 7i ^o 

< ^iiN £(<^) J+1 ^ «iM- 

□ 

We now turn the attention to the nonlinear term. 

Proposition 4.5. Let E G OPC7- m '- n (M d ), m > 1, n > 0, and h e N d , \h\ < 
max{n - 1, 0}, p fe G N d ; < m - 1, for 1 < k < I, I > 2. 

Then for every s > d/2 + max k {\p k \} there exists a constant C s > such that, 
for every e small enough and u G 5(IR d ), the following estimates hold: 

0<|a| + |/3|<AT 11 h|M|J fc=l 

<C s 5(||n||i- 1 ^l 1 M + (5^ 1 M)') 

if n > 1 and 

(4 - i7) e ^tumiiw 1 ^ns''-)!- 

< c^di^^Tiiiir^iiM + (££>])')■ 

i/0 < n< 1. 

Proof. Let n > 1, (hence < n — 1). In the sum (14.161) we consider the terms 
with a = 0. Namely, we prove that 

(4.i8) Yi M ii^(^n^)ii s <^ii«ii^ £ -iM- 

Given /3 ^ 0, let j G {1, . . . , d} such that ^ 0. Since £o2c i 2C h G OPGf" m '- r,+ l h l +1 (R d ) 
is bounded H s ~ M (R d ) -> # s (R d ), with M = max fc {|p fc |} (because < n - 1, 
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\pk\ < m ~ 1); an d applying Schauder's estimates (recall that s — M > d/2) we 
have 

_||£;(^n^^)ll s < ^ £ ( |^|_i)| ll x ^^IUmHIt 1 - 

Then (14.181) follows by writing 

(4.19) x^d^u= Y (- 1 ) l7l (/? ~ e,)l g pi _ 7( ^_ e ._ 7 ^ 

7<Pl 

in view of (I2.5p . and using 

(4 20) < ! 

1 • j (|/3|-l)!(/9- ei - T )! " (|/3|-1- | T |)!' 

cf. (J23D- 

For < n < 1, consider first the terms in (I4.17P with a = 0. We prove that 

0^|/3|<V fc=l 

To this end, given (3 ^ 0, let j G {1, . . . , d} such that /3j ^ 0. Since is bounded 
H s - M (R d ) -> # s (R d ), M = max fc {|p fc |}, by Schauder's estimates we have 



cf. the action on weighted Sobolev spaces described in Section [2j Then the claim 
follows by applying again (I4.19P and (I4.20p . 

We now treat the terms with a ^ in the sums (I4.16P and (14.171) . Namely, we 
prove that (both in the cases < n < 1 and n > 1) 

(4.21) £ m axlbl m ^ X da ^ X U dPku ))\U < ^(S^M) 1 . 

a<\ a \ + \f>\<N md *tl«b IPIJ-- k=l 

Let q^O and j e {1, . . . , d} such that ay 7^ 0. We can write 



^d a (x h Y[d" k u) = Q^[u]+Q^[u], 



x 

k=l 
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with 



Q? p [u] = d Xj x?&*- e > (x h \[d pk u), 



k=i 



Q^[u] = -fyxP-^d"-^ (x h Yl d pk u) . 



k=i 



Now we estimate ^~^~j^|yi ll-^Qi'^M lis- To this end observe that, by Leibniz' 
formula, 

S +Si + ... + Si = a-ej v ' k=l 

5 <h 

Let now 5 be a multi-index of maximal length among those satisfying | <5 1 < |^o| 
and 5 < (3. Observe that Ed Xja o x h - s °x' s ° e OPG~ m+1 >-»+l ft l(R d ) is bounded 
H s - M {R d ) ->■ H s {R d ) with M = max fe {|p fe |} (because |p fc | < m - 1 for 1 < fc < I 
and < max{n — 1, 0} < n). Hence 

>l + l/3| 

rll^MII- 



max{|a|, \/3\}\ 



SQ+6i + ... + Si=a-ej 
6 () <h 



£ N+W\ (a- ej )\ h\ 
max{|a|,|^|}!5 !5i!...5/!(/i-5o)! 



x \\x p ~ So Wd 5k+pk u\\ s - M - 



k=l 

We can now write 

i i 

(4.22) x p -~ s ° \\ d 5k+pk u = Yl x Jk d 5k+pk u, 

k=l k=l 

where 7i + . . . + 7; = /3 — <5o and, if \/3\ < \a\ — 1, with |7^| < \Sk\ for 1 < k < I (which 
is possible because in that case \/3 — <5o| < |« — <5o | — 1; observe that if \8q\ < \5q\ 
then (3 — 5 = 0), whereas, if \(3\ > \a\, with |7 fe | > \5k\ for 1 < k < I (which is 
possible because in that case \Sq\ = \Sq\ and \j3 — Sq\ > \a — Sq\ > \a — Sq\ — 1). 
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Hence we get by Schauder's estimates 



M+\P\ 

(4.23) ri , lon( \\EQf[u 



max{|a;|, \/3\} 



1 Aik\+\Sk\ 



max{ 1 7fc 1 , 1 5h I i ! " 

Jo+iX+.-.+^a-e^ fe=l Ll ,/C| ' 1 K|J 



where if |/3| < |a| — 1 we used the inequality 

( ' } {\a\-\)\6 Q \6 1 \...6 l \ ~ |5 |!|5i|!...|5i|!' 

which is (12.31) . whereas if \a\ < we applied 

(4-25) 1 < ^ 1 



|/3|U„!...<y - |J |!| 7l |!...| 7 ,|!' 

which also follows at once from (I2.3p . 

Now, we write x lk d &k+pk u = d Pk (x lk d Sk VL) + [x" /k d 5k , d pk ]u in the last term of 
ffl~23j) . so that 

|| z i*9**+« u || s _ M < \\x lk d 6k u\\ s + \\(D)- lpkl [x lk d 5k ,d pk }u\\ s . 
Using this last estimate we get 

rll^MII- 



max{|a|, 

,=-[7*1+1**1 



£ n j { | 7t |,| it | } i {«^ e ik^i^.^ii.}, 

S Q+ 4 1+ ...+i J= a- ej fc=l «IJ |7|<m-l 

<5 <'i 

(recall that the 7fc's depend on (3,5i, . . . ,5{). We now sum the above expression 
over | a | + \/3\ < N, a =^ 0. When a and /3 vary, every term in the above product 
also appears in the development of 

E majl{|5| E iiw-ww.} }, 

|&|+|£|<JV-1 mdX tl«l' IPI/- |7l<m-l 
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and is repeated at most, say, L times, with L depending only on h and the dimension 
d. Hence we obtain 

E 7i i \m EQ " Nil. 

|S|+|£|<jV-l max ll a l' l^l/- M<m-1 

< C':e{S s /_M + C':'eS s /_ 2 [u]} 1 < PT^iM)'. 

where we used Proposition 14.41 applied with <9 7 and (D) - ' 7 ' in place of P and E 
respectively, and we understand S'lffu] = 0. 
We finally show that 

F \a\+\P\ 

£ ^7^II^MII.<W^H)'. 

\a\ + \P\<N U " 1 " I J 

Since the arguments are similar to the previous ones, we give only a sketch of the 
proof. 

We can write 



n a 'P\ u ] -6- V (a-ej)\ h\ j-^tt. 

« +«i+...+ff I =a-e J - U 1 ( V U7 fc=1 

<5 <h 



If /3j 7^ 0, we choose a multi- index 5q of maximal length among those satisfy- 
ing \5 Q \ < \5 Q \ and 5 < (3 — ej. Next, we use the fact that E o x h ~ 5 °x s ° e 
0PG - m ,- n+ |ft|( R ^ ig bounded with M = max fc {|pfc|} and 

we apply the decomposition 

i i 

k=l k=l 

with 7! + . . . + 7; = (3 — ej — So, and | o^fc I < |^| for 1 < k < I, if \(3\ < \a\, or 
\lk\ > l^fcl for 1 < k < I, if \a\ < \(3\. Moreover, if |/3| < \a\ one uses 

Pi 1 (Q~e 3 -)! 1 
|a|(|a|-l)!<J Wi!...W " WW . . . \5 t \\ ' 

in place of (14. 24 p . whereas if \a\ > \(3\ (hence \S \ = \S \) one uses 
j 1 ( a - ei )! < _ l 



|/9|O0|-1)! do!...*,! " |5 |!|7i|!---b|! 
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in place of (14.251) . Therefore we get the same formula (14.231) with in place of 
Qi- The proof then proceeds as that for Q± without other modifications. □ 

We are now ready to conclude the proof of Theorem 14.11 

End of the proof of Theorem \4-l\ (the case m > 1). It follows from ( 14. ip that, for 
a, (3 e N d , e > 0, 

e M+lfl £ N+W\ e \«\+W 

X ^d a Pu = r , | lon A / + n i t on, x QaF H 



max{|a|, max{|a|, \/3\}\ max{|a|, 
so that 

e i«i+ii8i £ \c\+m £ w+\p\ 

max{|a|, |p|}! max{|a|, |p|}! max{|a|, |p|}! 

max{|a|, |p|}! 

We now apply to both sides the parametrix E of P. With i? = EP — I G 
OPG- 1 - 1 (M d ) we get 



x f } d a u = - T ---—R(x f3 d a u) + - T --——E[P,x p d a }u 



max{|o;|, \/3\}\ max{|a|, \/3\}\ max{|a|, 



£ W+\P\ £ \<*\+\P\ 
max{|a|, \p\\\ max{|a|, \p\\l 

Taking the H s norms and summing over \a\ + \/3\ < N give 

(4.27) ffl < IMI.+ E S^Sffli W^II- 

0<|a|+| / 8|<JV 

+ y — £ n"r la nM E (^d a f)i 

0<M + |/3|<V 
0<M + |/3|<Af 

The second and the third term in the right-hand side of (14 . 2 7[) can be estimated 
using Propositions 14.31 and 14.41 while the term containing / is obviously dominated 
by S%£[f]. For the last term we can apply Proposition 14.51 Hence, for n > 1, we 



max{ a 


M\Y 


e [«l+ 


\p\ 


max{ a 


M\Y 




m 


max{ a 


M\Y- 
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have that, for e small enough, 

s s /[u] < \\u\\ s + c s s^if] + c#(£%lm + (O^M) 1 + IMI^LiM)) , 

whereas if < n < 1 we get 

S s /[u] < \\u\\ s + C.S£[/] + C s e(s s /_M + £ ((^liN)' 

i 

+ \\(x)^u\\ l ^S s I f_ 1 [u])). 

In both cases we obtain 5^f[w] < oo if e is small enough, which implies u G 
AsectQ& d ) by Proposition[33](or, more simply, by the standard Sobolev embeddings, 
since s > d/2). 

□ 



4.2. Proof of Theorem I4.lt the case < m < 1. In this case the nonlinearity 
has the form 

(4.28) F[u] = Y,F h ,ix h u\ 

h,l 

where I G N, I > 2, h G N d , with < maxjn — 1,0} and i*}^ G C, the above sum 
being finite. 

We follow the same argument used for the case m > 1. In particular we can 
estimate the first four terms in the right-hand side of (14.271) as before since Propo- 
sitions H]3] and S3] hold in general for m > 0. Hence, to conclude, it is sufficient to 
prove an estimate for the nonlinear term. We have the following result. 

Proposition 4.6. Let P satisfy the assumptions of Theorem \4-l for < m < 1 

and let E be a parametrix of P. Let I G N, / > 2, h G N d , \h\ < max{n — 1,0}. 
Then, there exists a constant C' s > and, for every r > 0, there exists C T > such 
that, for every e small enough and u G 5(R d ) we have 



(4.29) + C',{eC T + t + + C',4{x) — u\\'-'S 



l-l c s,e 
N-l 



Proof. We first consider the terms in (I4.29[) with a = 0. Since Eox h G OPG~ m '°(R d ) 
is bounded on H s {R d ), we have, by Schauder's estimates: 



\E(x h+0 u l )\\ s < C' a \\xPu l \\ B < C'^\\x^u\\ s • \\ Xj u 
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if, say, (3j ^ 0. Then we get 



(4.30) mII^^II^^IK^^II^-^-iM 

0<\P\<N 

Consider now the terms in (14. 291) with a^O. We may write 



Pd a {x h u l ) = x h+p d a (u l )+ ( h ) ( a ^ v x h+ ^d a -\u l ) 
= x h d j {xPd a ~ e * (u 1 )) - p j x h+f> - e *& x - e * (u l ) 
7/ (a-7) ! * 



(4.3i) + z {':)t^x,* m -^-^')- 



Since E is bounded H s - m (R d ) H'W d ), we then obtain 



\\E(x^d a (x h u l ))\\ s < C'^djix^d^iu^Ws^ + C'JjWx^d^iu 1 ) 

7/ (a -7) ! 



(4.32) + c: £ (^J^y^^u- 



j<h 

Let us estimate the first term in the right-hand side of (I4.32p . We observe that for 
every r > there exists a constant C T > such that 

Hence 

11^(^9°-^ (U^ll-m = ||(£>)- m ^(x^-^(^))|| S 

< r||a,-(a^a°'- e '(« , ))|| s + C7 T ||^9«-^(« , )||. 

< r / 8 i ||a^- e '9 a -^(u , )||. + r||x /3 9 a (u / ) ||, 
+C T ||^9 a -^(«')|| s . 

Now we replace 

! 1 

d a {u l ) = lu l ~ l d a u + V — - - TT d Sk u 

«!+... +5,=a 1 1 fc=l 
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in the last estimate and we come back to ( 14.321) . We get, for a new constant C' s > 0, 

\\E(x^d a (x h u l ))\\ s < C:C T \\x^d a -^(u l )\\ s + C:(l + r)/3 3 \\x^ e W a -^(u l )\U 

+TC'\\u\\ l - 1 \\x^d a u\\ s 



+< E ra'^n^i 

d]_ + ... + d^= CK 

& k ^a Vfc 



(4.33) +C ; £ g)^L^a~V 



7<A 



We have now to estimate the terms in the right-hand side of f)4.33p . Concerning 
the first one, applying Leibniz' formula we obtain 



6 1 \...6 l \ 

oi+...+<5;=a— e-j k=l 

If |/3 1 > \a\, then we can argue as in the previous section and find ji, ■ ■ ■ ,ji G N d 
such that 71 + . . . + ji = ft and \jk\ > \6k\ for every k = 1, . . . , I. Moreover, we 
observe that the following estimate holds: 

1 (a — eA\ i 
( 4 -34) ToTT • T~i TT < 



\ft\\ 6 1 \...6 l \ ~ | 7 i|!...|7i|l' 
Then 
(4.35) 

_ M+\P\ _ _ ±,. £ \ik\Mh\ 

e ^w7w l|i S °" ei(u)ll '- c " £ £ £ n— 

\<*\ + \f)\<N LI DI^IJ |a| + |/3|<JV 5i + ...+,5 ; =a-e, fc=l 1 ,B| 

0<|a|<|/3] 0<|a|<|j8| 

x \\x^d 5 <°u\\ s < G": e(S S rf_MY- 

On the other hand, for \ft\ < \oc\ — 1 we can choose multi-indices ji, ... ,ji such 
that 7i + . . . + ji = ft and |7fc| < \6k\ for any k = 1, . . . ,1 and observe that 

(4-36, < 1 



(H -!)!<?!!...<?,! " \6t\ 
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Then 

_ F \o\+m 

< 4 ' 37 > E „„/!„, Ifl n, ll^°-' J ("')II, 



max-flal, 1/31}! 

0<|/3|<|a|-l 



' £ hk\+\Sk\ 



< c " E nno~ l|x ^ tt|l *-^ /£(,s ^ i[ul)J ' 

| a |+|/3|<jv 5i+...+5i=Qi-ew fc=l ' 
0<|/9|<|q|-1 

for new constants C" and C"'. 

For the second term in the right-hand side of (14 . 3 3 j) we can argue as before, with 
7i + . . . + ji = (3 - ej and |7 fc | < \5 k \ for k = 1, if |/3| < \a\ or |7 fc | > \5 k \ for 
k — 1, . . . , I, if |/3| > | or | , using the estimates 

& 1 (a-e i )! < 1 



1/31(1/31-1)!^!...^! " | 7 i|!...|7/|! 
respectively, 

Pi 1 (a-e 3 -)! 1 
|a| (|a|-l)! 5i!...5«! ~ |5i|!...|^|! 
instead of (I4.34p . respectively (14.361) . We obtain, for a new constant C' s > 0, 
^ £ \<*\+\P\ 

<« 8 > H E s =snOT»^^^»- £ 

Concerning the fourth term in (I4.33P , we can decompose similarly /3 = 71 + . . . + 7; 
and argue as before, taking into account that now I71+. . .+ji+5i + . . .+5i\ = \a+(3\, 
so that we do not longer gain e as a factor in the estimate. Hence we get 

{4M > e max{ | ttU/3 | }! e ^i^n^ii.s^^M)' 

|a+/3|<JV u !' 1 ' 1 J «!+•■.+«/=« fc=l 

Q^O 4 fc ^QfVfe 

for a new constant C^. Finally, for the last term in (14 . 3 3 [) . we first observe that 
max{|a — 7], |/3 — 7]} = max{|a|, |/3|} — |7|. Then we can argue as before obtaining 
the estimate 
(4.40) 

i.^iv max{|a|, |/3|}! ^ (« - 7)! 

7 </j, 

The estimates (00) . fl435|) . (Q7J) . pg| . (fQ5j) . fl440|) applied in pg| yield 



HOLOMORPHIC EXTENSION FOR ELLIPTIC EQUATIONS 27 

End of the proof of Theorem \4-l\ ( the case < m < 1). Using the same argument 
as in the case m > 1, by Propositions 13.31 14.3} I4.4[ 14.61 we obtain 

S s /[u] < \\u\\ s + C' s S^[f) + C^tM + (rC'MVS'/iu} 

i 

+ C' s (sC T + r + eXSj^M)' + C'^Ux^ut-'S^M) 

for every N > 1 and e small enough. Now, choosing r < (2 Csll^Hi -1 )~ 1 we 
obtain 

S^>] < 2||«||. + 2C' s S s J[f] + 2C' s eS s /_ 1 [u]+ 

+ (2C' s (eC T + r + £)(S s /_ 1 [u]) 1 + 2Cie\\(x)^ u \fc 1 S s /_M)- 
i 

Then we can iterate the last estimate observing that, shrinking r and then s, the 
quantity eC T + r + e can be taken arbitrarily small. This gives Sjf [w] < oo and 
therefore u G A sec t{^ d )- 

5. Remarks and applications 

5.1. Lower a priori regularity. For special nonlinearities, in Theorem [LT] we can 
assume lower a priori regularity on the solution u. For example, if F[u] = (d p u) 1 , 
\p\ < min{m — 1, 0}, / G N, / > 2, or even F[u] = \d p u\ l ~ 1 d p u, \p\ < min{m — 1, 0}, 
I G N, I > 2 odd (as in (02) ). then we can assume u G H 3 (R d ), with s > | - 2z^t, 
s > Indeed, such a solution is actually in if 00 (R a! ); see e.g. [5j Lemma 4.1, 
Remark 4.1] (where that threshold is also proved to be sharp). We also refer to [5] 
for other types of non-linearities. 

5.2. Eigenfunctions of G-elliptic operators. In the linear case, the assump- 
tions on the a priori regularity of u can be relaxed assuming u G S'(M. d ). Then, we 
have the following result. 

Theorem 5.1. Let P be a G- elliptic pseudodifferential operator with a symbol 
p{x,^) satisfying (11.91) . Then there exists e > such that every solution u G <S'(R d ) 
of the equation Pu = extends to a holomorphic function in the sector of C d 

C £ = {z = x + iy G C d : \y\ < s(l + \x\)}, 

satisfying there the estimates (jl.3p for some constants C > 0, c > 0. 

Proof. It follows from the existence of a parametrix (see Section [2]) that any so- 
lution u G iS'(M d ) of Pu = is in fact a Schwartz function. Hence we can apply 
Theorem 11.11 directly without using Lemma 14.21 Moreover it follows from the clas- 
sical Fredholm theory of globally regular operators (see e.g. [271 Theorem 3.1.6]) 
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that the kernel of P is a finite dimensional subspace of 5(IR d ), which implies that 
there exists a sector where all the solutions extend holomorphically. □ 

The main application of Theorem 15.11 concerns eigenfunctions of G-elliptic oper- 
ators of orders m > 0, n > 0. Indeed, in that case, if P is G-elliptic also P — A 
is G-elliptic for every A G C, and one can apply Theorem 15.11 to P — A. As re- 
gards existence, we recall that if P G OPG m,n (R d ), m > 0, n > 0, is formally 
self-adjoint (i.e. symmetric when regarded as an operator in L 2 (M. d ) with domain 
<S(IR d )) then it has a sequence of real eigenvalues either diverging to +oo or — oo, 
and L 2 (M d ) has an orthonormal basis made of eigenfunctions of P (cf. e.g [25] or 
pTf Theorem 4.2.9]). As an example in dimension 1, one can consider the operator 
Pu = —(1 + x 2 )u + x 2 u — 2xu', i6R. 

5.3. Solitary waves. The present subsection is devoted to some applications to 
solitary waves, in particular to the proof of Theorem 11.21 First we report the 
following useful characterization of the condition fll . 16j) . 

Proposition 5.2. The estimates f l 1 . 1 6 [) are equivalent to requiring thatp(^) extends 
to a holomorphic function p(^ + ir]) in a sector of the type (I1.18p . and satisfies there 
the bound \p(£ + irj)\ < C"(£) m . 

Proof. The sufficiency of f ll . 16[) for the holomorphic extension with the desired 
bound follows exactly as in the last part of the proof of Theorem I3.2| where in 
(13. 5 p the exponential factor is now replaced by (£) m . 

In the opposite direction, we obtain f 1 1 . 1 6 j) from Cauchy's estimates applied to a 
disc in C with center at £ and radius for some small e' > (independent of 

0- ' □ 

Proof of theorem Consider first the equation (11.141) . We observe that u satisfies 
the equation 

(5.1) Mu + (V - l)u = F[u], 

cf. the proof of [91 Theorem 3.2.1]. By (I1.17P and the condition V > 1 the symbol 
of the linear part of the equation (15. ip . that is + V — 1, is G-elliptic: for some 
constant c > 

(5.2) p(0 + V - 1 > c{0 m , (6t 

Moreover by f 1 1 . 1 6 [) it satisfies the analytic symbol estimates f ll .Qp (with n — 0). To 
conclude the proof it is sufficient to show that u G H s (M. d ) and (x) £ °u G L 2 (M. d ) for 
some eo > 0, s > d/2, and to apply Theorem ll.il The fact that u enjoys the above 
properties will follow from [91 Theorem 3.1.2, Corollary 4.1.6] once we observe that 
the function 

m) 1 



p(0 + v-i 
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satisfies \K(g)\ < C{£)- m for some C > and belongs to H°°{R d ). This is clear, 
because (15. 2p and ( ll.lfip give 



1^(01 <c«<0 



-m— |a| 



and m > 1. 

The case of the equation ( 11.15)) is completely similar: in place of (15. 2p one 
just has V Mu + V — 1 = F[u], and the above arguments apply to the function 

m) = (Vp(o + v-iy 1 . 

The theorem is then proved. □ 

As an example where the solutions are known in closed form, consider the gen- 
eralized Korteweg-de Vries equation 

(5.3) v t + v x + v l v x + v xxx = 0, 

where / > 1 is a positive integer. Here we have = £ 2 . The solitary wave 
solutions have the form v(x, t) = u(x — Vt), where V > 1 and 



which has poles at the points z = ijj^=f, k £ Z. Also, the exponential decay in 
sectors containing the real axis predicted by Theorem 11.21 is confirmed. 

During the years 1990-2000, several papers were devoted to 5-th order and 7- 
th order generalization of KdV, see for example Porubov [291 Chapter 1]. The 
corresponding stationary equation is of the type 

m 

(5.4) ^2cLjU {j) + Q[u] = 0, 

j=0 

where Q is a polynomial, Q[u] = J2jL2^j u ' j anc ^ a o 0- Because of physical as- 
sumptions, the equation J2 1 jLo a i^ = ® nas no P urer y imaginary roots, and then all 
the solutions of the corresponding linear equation have exponential decay/growth. 
This condition can be read as G-ellipticity of the symbol of the linear part of 
the corresponding stationary equation: ^JL a i(^)' ? 7^ for £ G M, in particular 
£ 2 + V — 1 7^ in the case of (15. 3p . Non-trivial solutions u of (15 ,4p with u(x) — > as 
x — > ±oo may exist or not, according to the coefficients dj, bj, and when they exist, 
in general they do not have an explicit analytic expression. Holomorphic extension 
and exponential decay on a sector are granted anyhow by Theorem 11.21 

5.4. Standing wave solutions of the Schrodinger equation. Consider the 
Schrodinger equation in 

id t v + Av = /iM*~ V (t, x)eRy Tn "' 
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with I G N, I > 2 odd, /i 6 C, and look at standing wave solutions, i.e. v(t, x) = 
e luJt u(x), u > 0. The corresponding equation for u is 

Au — cou = fi\u\ u. 

Since the operator A — u is G-elliptic (because cu > 0), when solutions u exist, with 
u G F s (M d ), s > | - s > 0, and (x) e °w G L 2 (IR d ), for some e > 0, then Theo- 
rem 11.11 and the remark in Subsection 15.11 assure that u extends to a holomorphic 
function on a sector of the type (11.131) and displays there an exponential decay of 
type (jl.3p . This applies, in particular, to the bound states in H 1 {R d ) exhibited in 
[I] when I < g§ , d > 3. 

5.5. Sharpness of the results. Here we show the sharpness of Theorem 11.11 as 
far as the shape of the domain of holomorphic extension is concerned. 
Consider in dimension d — 1 the equation 

— u" + e- 2ie u = ^-n 2 , 

where — n < 9 < n, \9\ ^ |. This equation is G-elliptic, since it is elliptic and 
£ 2 + e~ 2ie 7^ for every (el. An explicit Schwartz solution is given by 

u{x) = 3Cosh~ 2 y— — %J ■ 

The function u extends to a meromorphic function in the complex plane with poles 
at z = e^ £l+7r//2 - ) (2k + 1)tc, k G Z. This shows that in Theorem 1 1.1 1 we cannot replace 
the sector f ll,13p . e.g., with a larger set of the type 

{z = x + iyeC d : \y\ < e(l + \x\)ip(x)}, 

for any continuous function ijj(x) > 0, with if)(x) — > +oo as \x\ —> +oo. 
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